X(3872) — > J/?pn + n~~ AND X(3872) — > J/?pit + n 7T° DECAY WIDTHS FROM QCD SUM RULES 



F.S. Navarra and M. Nielsen 
Instituto de Fisica, Universidade de Sao Paulo, CP. 66318, 05389-970 Sao Paulo, SP, Brazil 

New spectroscopy from the B factories, the advent of CLEO-c and the BES upgrade renewed 
the interest in charmonia. Among the new measurements, the state X(3872) has received special 
attention due to its unexpected properties. Its structure has been studied with different theoretical 
approaches, most of them being able to reproduce the measured mass. A further test for the 
theoretical descriptions of the X(3872) is to explain its narrow decay width. In this work we address 
the decays X — > J/ifj 7r + 7r _ 7r° and X — > J/ip -k + it~ , using QCD sum rules with the hypothesis that 
A is a four quark state. 



PACS numbers: 11.55.Hx, 12.38.Lg , 12.39.-x 



During the last three years several new hadronic states have been observed, as for example the D^j(2317) pj and 
the A (3872) 0. The experimental observations were always followed by theoretical efforts to understand the basic 
properties of the new particles, in particular the mass and the decay width. In the charm sector, simple potential 



models, which had been so successful in the past, failed in reproducing the masses of the new states. This was taken 
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as an indication that these particles are not simple q — q bound states. As for the very narrow decay width, whereas 
it is expected in the case of the £)j^(2317), since it decays through an isospin violating channel, in the case of the 
A (3872) it is really surprising. Even more surprising is the observation, reported by the BELLE collaboration 



that the X decays to J/tp 7r + 7r~7r°, with a strength that is compatible to that of the J/ipTr + 7r~ mode: 

Br(X -► JMtt+tt-tt ) 

■ „\ -n—— r = 1.0 ±0.4 ±0.3 . (1 

m ■ 

This decay suggests an appreciable transition rate to J/tpu; and establishes strong isospin violating effects. The 
measured A (3872) mass can be reproduced in several approaches and it is not yet possible to discriminate between 
the different structures proposed for this state: tetraquark 0,O|, cusp ID, hybrid 0, glueball y| or DD* bound state 

*■■»«. ' The theoretical study of the decay width can help in clarifying this situation. In this work we use the method 
^ i of QCD sum rules (QCDSR) III El to studythe hadronic decays of X(3872) given in Eq.©, considering X 
^f^- as a four-quark state. In recent calculations 0,0,0, the QCDSR approach was used to study the light scalar 
Q_i| mesons, the -D^j(2317) meson and the X(3872) meson considered as four-quark states and a good agreement with the 
experimental masses was obtained. In particular, in ref.^3 we have considered the A(3872) as the J PC = 1 ++ state 
'-^ ' with the symmetric spin distribution: [cq\s=i[cq]s=o + [cq]s=o[cq]s=i- The interpolating field for X q is given by: 



i« = !!^£[(gJC 75 c fe )fe 7 ,CcD + (q^C lfl c b )(q dl5 CcJ)] , (2) 

where a, 6, c, ... are colour indices, C is the charge conjugation matrix and q represents the quark u or d. 

As pointed out in isospin forbidden decays are possible if X is not a pure isospin state. Pure isospin states are: 

X(/ = 0) = ^-^, and X(I=1) = ^J^. (3) 

If the physical states are just X u or X4, the mass eigenstates, maximal isospin violations are possible. Deviations 
from these two ideal situations are described by a mixing angle between X u and Xd 0: 

Xi = X u cos 6 + Xds'md, 

X h = -X u sm6 + X d cos9. (4) 

In ref. 4] , by considering the ratio of branching ratios given in Eq.Q, they arrived at 9 ~ 20° and at T(X — > 
J/^7T7r) ~ 5 MeV. However, to arrive at such small decay width they had to make a bold guess about the order 
of magnitude of the XJ/ipV (where V stands for the p or ui vector meson) coupling constant: gx^pv — 0.475. In 
this work we evaluate the XJ/ipV coupling constant directly from the QCD sum rules. For the light scalar mesons, 
considered as diquark-antidiquark states, the study of their vertex functions using the QCDSR approach was done 
in ref. 17]. The hadronic couplings determined in ref.^3 are consistent with existing experimental data. In the case 
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of the meson (2317) considered as a four-quark state, the QCDSR evaluation of the hadronic coupling constant 
9d sJ d s -k HI gives a partial decay width in the range 0.2 keV < r(£>+ 7 (2317) -> D+n°) < 40keV. 

The QCDSR calculation for the vertex, X(3872) J/ipV, centers around the three-point function given by 

rWp,p',<7) = J d 4 xd 4 y e** e^yn^ a (x,y), with U^ a (x,y) = {0\T\fi(x)fi{y)jZ\o)]\0), (5) 

where p = p' + q and the interpolating fields arc given by: 

Of = CaJfiCa, (6) 



3v = ^{UalvUa + {-lYda^da), (7) 

with N = I = 1 for V = p, N = 1/3, I = for V = cj and 

£ = att + bj«, (8) 

where j£ is given in Eq.© and (see Eq.@). 

. „ | a = cos 6* r f a = — sin# , nS 

' i . f ° rJ M ( 9 ) 



& = sin ' [ b = cos 

Using the above definitions in Eq.JSJ) we arrive at 

— iN 

iWz >y ) = -^j= (a n^^.y) + (-1) 7 & n^ Q (x,y)) , (10) 

with 

n^Jz, y) = ea6 C e dec Tr [5 e c a , b (x) 7p ^, b ( a :)( 7 5C^(y)C7,C^(-2/)C 7ct - 7a C^(y)C 7 ,C^(-y)C 7 5)] , (11) 

where S% b (x — y) = (0\T[q a (x)qt(y)]\0) is the full quark q propagator. 

To evaluate the phenomenological side of the sum rule we insert, in Eq.JSJ, intermediate states for X, J/ip and V . 
Using the definitions: 

(0|$| J/rPip')) = m+Uertf), (0\tf\V(q)) = m v f v e u (q), (X q (p)\f a \0) = \ q e* a (p), (12) 
we obtain the following relation: 

(p 2 -m2 f )(y-m,2)(g2_ m 2 7 ) ^ m ^ TO ^ 



where the dots stand for the contribution of all possible excited states, and the form factor, gxipvio 2 ), is defined by 
the generalization of the on-mass-shell matrix element, (J/ifjV\X), for an off-shell V meson: 

which can be extracted from the effective Lagrangian that describes the coupling between two vector mesons and one 
axial vector meson Q : 

C = %9xi > v^ vaa ^X v )^ a V a . (15) 

From Eq. (|13f) we see that we have four independent structures in the phenomenological side. For each one of these 
structures, i, we can write 

iq ( P hen) Aj gx^vjq 2 ) f°° P C i 0nt (p 2 ,q 2 ,u) 



(q 2 - niy)(p 2 - m 2 x ){p' 2 - mV) J± m * u - p' 2 
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In Eg. (|16f) . pf nt {p 2 ,q 2 ,u), gives the continuum contributions, which can be parametrized as pf nt (p 2 , q 2 , u) = 
b% J^lp^ Q( M ~ u o) l2ll with so and u$ being the continuum thresholds for X and J/ip respectively. Tak- 
ing the limit p 2 = p' 2 = —P 2 and performing a single Borel transformation to P 2 — ► M 2 , we get (Q 2 = —q 2 ): 

nf^(M 2 ) = ^sxMQ 2 ) L- m %/M> _ e -m° x/ M>\ b e _ So/M , r~ e -u/M> d 

(m 2 v + Q 2 )(m 2 x - m\) V J (m\ + Q 2 ) Ju a 

.(17) 

where P>i and pl c (u,Q 2 ) stand for the pole-continuum transitions and pure continuum contributions. For simplicity, 
one assumes that the pure continuum contribution to the spectral density, p1 c (u, Q 2 ), is given by the result obtained 
in the OPE side for the structure i. Asymptotic freedom ensures this equivalence for sufficiently large u. Therefore, 
one uses the Ansatz: p1 c {u,Q 2 ) = pf PE (u,Q 2 ). In Eq.l|17|). Bi is a parameter which, together with the form factor, 
gxipv(Q 2 ), has to be determined from the sum rule. 




(c) (d) (e) 



FIG. 1: Diagrams which contribute to the OPE side of the sum rule. 



In the OPE side we work at leading order and consider the condensates up to dimension five, as shown in Fig. 1. 
To keep the charm quark mass finite, we use the momentum-space expression for the charm quark propagator. We 
calculate the light quark part of the correlation function in the coordinate-space, which is then Fourier transformed to 
the momentum space in D dimensions. The resulting light-quark part is combined with the charm-quark part before 
it is dimensionally regularized at D = 4. For each structure i, we can write the Borel transform of the correlation 
function in the OPE side in terms of a dispersion relation: 



n: 



q(OPE) 



(M 2 ,Q 2 ) 



4 m' 2 



A 



q{OPE) 



(u,Q 2 ) 



-u/M' 



du 



(18) 



where the spectral density, p q ^ OPE ^ , is given by the imaginary part of the correlation function. The perturbative term 
(diagram in Fig. 1(a)) contributes only to the structures e «/ iI/ < 7 p^ and t aiLai p' a q 1 q Vl while the quark condensate and 
mixed condensate (diagrams (b) to (e) in Fig. 1) contribute to the structures t a ^ va q a and e Qt/cr7 p^g 7 p^. Therefore, 
to get more terms contributing in the OPE side we have two options for the structures: e a ^ va q IJ and e a " a " r p' (J q 1 p^. In 
order to test the dependence of the results with the chosen structure, we will work with these two structures. 

Transferring the pure continuum contribution to the OPE side we get for the structure e a<y<7 " 1 'p'^q^p' (which we call 
structure 1): 



Hi 



q(OPE) 



(M 2 ,Q 2 ) 



■i(qq) 

3tt 2 Q 2 



to5 
2 5 



4m? 



da 



l + 3a 



e o(l-a)Af^ 



(19) 



and for the structure t a ^ va q a (which we call structure 2) we get: 



uf OPE) {M 2 ,Q 2 



i{qq) 

3tt 2 Q 2 



3Q 2 



du e~ u ' M uJl-Am 2 /u - + -± | + 

1 2 ? 



4 



2 5 



j l + 3a 
da ~2n \ 



(20) 



In Eqs. I|19|) and (|20[1 we have used the relation (qga.Gq) = m^qq). 

Making use of Eqs. 10 and i|l(J|) . and working at the SU(2) limit, i.e., considering the quarks u and d degenerate, 
we arrive at three sum rules for each structure, that can be written in the general expression: 



where 



and 



C 



xv 



(Q 2 ) (< 



/M2 )+B t e- s «/ M2 



c? v (Q 2 



-m%/M 2 



^ 2 K 2 A XV (Q 2 ), and C? v (Q 2 
m x - 



in 



-A XV {Q 2 



x 



(21) 



(22) 



Ax lP {Q 2 ) 
Ax hP {Q 2 ) 
AxUQ 2 ) 



cos 9 + sind 2 



COS f - 

cos( 



sin 
— sin 



gx^piQ ), 
9x h ^ P {Q 2 ), 



—rripfp- 

cos V + sin ( 



(23) 



Since from Eq. (|21|l we see that the OPE side of the sum rule determines only one value for C xv for each structure 
(for a fixed value of Q 2 ), we arrive at the following relations between the form factors: 



gxi^u{Q 2 ) 
gx^piQ 2 ) 
gx h ii>u>(Q 2 ) 
gx h i>p(Q 2 ) 



rripfp cos 9 + sin 9 
Smuifuj cos 9 — sin 9 ' 
Ttipj ' p cos 9 — sin 9 
3m w / w cos 9 + sin 9 



(24) 
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FIG. 2: Dots: the RHS of Eq.J5TJ, as a function of the Borel mass for Q 2 = 3 GeV 2 . The solid line gives the fit of the QCDSR 
results through the LHS of Eq.(|2TJ for As = 0.5 GeV. 



In the numerical analysis of the sum rules, the values used for the quark masses and condensates are: m c = 
1.2 GeV, (qq) = — (0.23) 3 GeV 3 , m 2 , = 0.8 GeV 2 . For the meson parameters we use their experimental values 
H3: m p = 0.776 GeV, m u = 0.782 GeV, m v , = 3.1 GeV, m x = 3.872 GeV, f p = 0.157 GeV, f u = 0.046 GeV and 
ftj, = 0.405 GeV. We evaluate our sum rules in the range 2.0 < M 2 < 3.0, which is the range where the two-point 
function for X(3872) shows good OPE convergence and where the pole contribution is bigger than the continuum 
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contribution |Tj|. We also use three different values for s = (3.872 ± As ) 2 GeV 2 : As = 0.4 GeV, As = 0.5 GeV 
and As = 0.6 GeV. For u Q we use u Q — (m^ + 0.5) 2 GeV 2 . The meson-current coupling, X q , defined in Ea. ((T2*|l . 
can be determined from the two-point sum rule [l^- I n Table I we give the results obtained from ref.0] for three 
different values of so. 

Table I: Numerical results for the meson-current coupling 



X q (GeV 5 ) 


Aso ( GeV) 


(1.85 ±0.01) x 10~ 2 


0.4 


(1.94 ±0.03) x 10~ 2 


0.5 


(2.02 ±0.06) x 10 -2 


0.6 



We start with the structure 1. In Fig. 2 we show, through the circles, the right-hand side (RHS) of Ea. l|21|l for 
Q 2 = 3 GeV 2 , as a function of the Borel mass. 

To determine gxipv(Q 2 ) we fit the QCDSR results with the analytical expression in the left-hand side (LHS) of 
Eq.(|23J, and we get (using As = 0.5 GeV): C? V (Q 2 = 3 GeV 2 ) = 7.01 x 10~ 4 GeV 5 and B x = -1.28 x f0~ 3 GeV 5 . 
Using the definition of Cf v (Q 2 ) in Eq.|22J) we get A XV {Q 2 = 3 GeV 2 ) = 1.49 GeV 2 . Allowing Q 2 to vary in the 
interval 2.5 < Q 2 < 4.5 GeV 2 , we show, in Fig. 3, through the circles, the momentum dependence of Axv(Q 2 )- 

From Ea. H23[) . we see that all form factors are related with the function Axv(Q 2 )- Since the coupling constant is 
defined as the value of the form factor at the meson pole: Q 2 — —my, to determine the coupling constant we have 
to extrapolate the QCDSR results to a Q 2 region where the sum rules are no longer valid (since the QCDSR results 
are valid in the deep Euclidian region). To do that we parametrize the QCDSR results through a analytical form. In 
Fig. 3 we also show that the Q 2 dependence of A X v{Q 2 ) can be well reproduced by the monopole parametrization 
(solid line): 

■ 4 -<« 2 ' = Q^rfn ■ < 25 » 

from where we can extract the value of AxviQ 2 ) & t the meson pole: AxviQ 2 = ~ m v) = 1-62 GeV 2 . 
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FIG. 3: Momentum dependence of Axv for Aso = 0.5 GeV. The solid line gives the parametrization of the QCDSR results 
(circles) through Eq. i2*5l . 

Doing the same kind of analysis for the other values of the continuum threshold we show, in Table II, the monople 
parametrizations of the QCDSR results, as well as their values at the off-shell meson pole. 

Table II: Monopole parametrization of the QCDSR results for the structure 1, for different values of Asq 



Aso (GeV) 


Axv(Q 2 ) (GeV 2 ) 


Axv{Q 2 = ~m 2 v ) (GeV 2 ) 


0.4 


70.2 
Q^+41.6 


1.71 


0.5 


66.8 
Q^+41.8 


1.62 


0.6 


63.8 
Q a +41.7 


1.55 
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FIG. 4: Dots: the RHS of Ea. l|21|l . for the structure 2, as a function of the Borel mass for Q 2 = 3 GeV 2 . The solid line gives 
the fit of the QCDSR results through the LHS of Eq.JH) for As = 0.4 GeV. 



In the case of the structure 2, the RHS of Ea. (|21() can also be very well parametrized with the analytical expression in 
the LHS of Eq.JHJ, as can be seen in Fig. 4. We get (using As = 0.4 GeV): C? v (Q 2 = 3 GeV 2 ) = 5.56 x 1CT 3 GeV 7 
and B 2 = -3.46 x 1(T 3 GeV 7 . Using the definition of C? v [Q 2 ) in Eq.® we get A X v{Q 2 = 3 GeV 2 ) = 1.29 GeV 2 . 
The Q 2 behaviour of Axv(Q 2 ) can also be well represented by a monopole form in the case of structure 2, with a 
precision similar to the one shown in Fig. 3. In Table III we give the monople parametrizations of the QCDSR results 
for the structure 2, as well as their values at the off-shell meson pole. 

Table III: Monopole parametrization of the QCDSR results for the structure 2, for different values of Asq 



As (GeV) 



Axv(Q 2 ) (GeV 2 



Axv{Q 2 



(GeV 2 



0.4 



59.0 
P^+42.6 



1.45 



0.5 



56.0 
Q^+42.: 



1.34 



0.6 



54.2 
3^+42.9 



1.28 



Comparing the results in Tables II and III we see that, althought the results from the structure 2 are somewhat 
smaller than the results from the structure 1, they are still compatible with each other. We will use these differences 
to estimate the uncertainties in our results. 

From Eq. (|23|l we see that, in the case of the meson cj, there is no mixing angle dependence in the relation between 
Axv an d gxipui- Therefore we can use the results in Tables II and III to directly estimate the XJ/ipuj coupling 
constant. We get 



9XipL 



= 13.8 ±2.0, 



(26) 

which is much bigger than the guess made in ref.0: gxipv = 0.475. 

Having the coupling constant and the relations in Eqs. I|23|l and l|24|) . we can estimate the decay widths of the 
processes X — » J / 'ip ir + tt~ tt° and X — > J/ifjir + ir~ by supposing that the 27r and 37r decays are dominated by the p 
and lo vector mesons respectively. In the narrow width approximation we have: 



dT 

ds 



(X -» J/1>(mr)) 



1 



%irm 2 x 



\M[' 



s Tymv 



p(s) 



2m 2 x 



(s - ml-) 2 + {m v T v f 



(27) 



with n = 2, 3 for V — p, U). In Eq. l|27|l . s is the invariant mass-squared of the pions, Ty and By-tnir are i respectively, 
the total decay width and the branching ratio of the V — > mr decay. The decay momentum p(s) is given by 



p(s) = 



y/A(m|,m 2 ,s) 



2m x 



(28) 



with A(o, 6, c) = a 2 +b 2 + c 2 - 2ab - 2ac - 26c. 



7 



The invariant amplitude squared can be obtained from the matrix element in Eo. l(14() . We get: 



\M\< 



2 

9xipv 



Am 



(rn 



,2 \2 



) , (' 



X 



2,s 



2ml 



(29) 



where we have replaced the form factor, gx-ipv(s) by the coupling constant gx-ipv, since from Tables II and III we 
can see that the form factor is very flat over the region, [n m„) 2 < s < (mx — to,/,) 2 , over which Ea. 1)27(1 will be 
integrated. Using the relations between the coupling constants from Eo. 1(24(1 we get the following relations between 
the decay widths: 



f T(Xi -> Jlj> 3tt) \ = m 2 p f 2 



\T{Xi -> J/V> 2tt) 
T(X h -> J/y, 3tt) 

rpo, -> j/v- 2tt) 



9mlfu; 2 



cosf 



sin! 



cosf 
cos f 



sm (V 
sin (9 



cosf 



smt 



V 



(30) 



where we have defined 
Iv = 



T v mv 



7T 



ds 



Am 



x 



+ 



2.s 



(m 2 x - s) 2 
2to 2 



p(s) 



2m^ 



(s - TOy ) 2 + (TOv-r y ) 5 



Since the relations in Eq. 1(301) do not depend on the value of the coupling constant we get 



fT(X Lh - J/V 3tt) 



0.152 



cos 9 ± sin 9 



\T(X Lh J/^i 2?r) / Vcosfy=Fsm( 
Therefore, using the central experimental data given in Eq.JQ) we obtain for the mixing angle 

9 ~ ±23.5° 



(31) 



(32) 



(33) 



for Xi or Xh respectively, which is in agreement with the result obtained in |4|: 9 ~ ±20°. 

Since, with the determination of the mixing angle in Ea. l(33|) by imposing the ratio in Eq.Q), we obtain the same 
width for any of the four decays in Eg. 1(301) . we can use the value of the coupling constant determined in Eg. ((26(1 to 
evaluate the partial decay width. We get 



T(X J/i> {nil)) = (50 ± 15) MeV, 



(34) 



which is much bigger than the experimental total width: T(X(3872)) < 2.3 MeV. 

As a matter of fact, a large partial decay width was expected in this case. The initial state already contains all the 
four quarks needed for the decay, and there is no violating rules prohibiting the decay. Therefore, the decay is allowed 
as in the case of the light scalars a and k studied in [17]], which widths are of the order of 400 MeV. However, even 
when there is no violatin g ru les prohibiting the decay, the decay can be prevented due to a non-trivial color structure 
in the initial state. In ref.[22|, an alternative technique was developed to obtain the form factor and coupling constant 
for multiquark particles. By multiquark we mean that the initial state contains the same number of valence quarks as 
the number of valence quarks in the final states. In this case, as can be seen in Fig. 1, the generic decay diagram in 
terms of quarks has two "petals", one associated with the J/ip and the other with the other vector meson V. Among 
the diagrams in Fig. 1 there are two distinct subsets. In the first (diagrams from (a) to (d)) there is no gluon line 
connecting the petals and, therefore, no color exchange between the two final mesons in the decay. A diagram of this 
type was called color-disconected (CD) diagrams in ref. [23. If there is no color exchange, the final state containing 
two color singlets was already present in the initial state. In this case the tetraquark had a component similar to 
a J/t/j — V molecule. The other subset of diagrams is represented by diagram in Fig. 1(e), where there is a color 
exchange between the petals. This type of diagram represents the case where the X is a genuine four-quark state with 
a complicated color structure. These diagrams are called color-conected (CC). In our approach we have considered 
all kinds of diagrams. However, if we consider only the CC diagrams, which means considering only the diagram (e) 

(CO 

in Fig. 1, we get g x ^ = 1-6 ± 0.3, and therefore 



CC 



(X -> J/ip (nir)) = (0.7 ± 0.2) MeV. 



(35) 
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This procedure may appear somewhat unjustified. However, we do believe that there should be a particular choice 
of the interpolating field, which represents a genuine four-quark state, for which CD diagrams vanish. From our 
calculation we find out that the interpolating field in Eq.Q has a component similar to a J/tp — V molecule. 

To summarize: we have presented a QCD sum rule study of the three-point functions of the hadronic decays of 
X(3872) meson, considered as a diquark antidiquark four quark state. Supposing that the physical state is a mixture 
between the isospin eigenstates, we find that the QCD sum rules result for the mixing angle is compatible with the 
result found in [4J. However, we get a partial decay width much bigger than the experimental total decay width. 
Therefore, we conclude that our particular choice of the interpolating field has a J/ipV molecule component, and is 
not the most appropriate candidate to explain the very small width of the meson AT (3872). Further studies, using 
different interpolating fields, are necessary for a better understanding of the structure of the meson A(3872). 
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